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I discuss challenge of lattice QCD, from quarks to nuclei, which connects QCD with nuclear physics. 



1. Introduction 

In 1935 Yukawa introduced a virtual parti- 
cle, now called the pion, to account for the nu- 
clear force , which bounds protons and neu- 
trons in nuclei. Since then the nucleon-nucleon 
(NN) interaction has extensively been investi- 
gated at low energies both theoretically and ex- 
perimentally. Fig. [TJ shows modern NN poten- 
tials, which are characterized by the following 
features |21 15] . At long distances (r > 2 fm ) 
there exists weak attraction, which is well un- 
derstood and is dominated by the one pion ex- 
change (OPE). At medium distances (1 fm < r < 
2 fm) , contributions from the exchange of multi- 
pions and/or heavy mesons such as p, to and a 
lead to slightly stronger attraction. At short dis- 
tances (r < 1 fm), on the other hand, attraction 
turns into repulsion, and it becomes stronger and 
stronger as r gets smaller, forming the strong re- 
pulsive core[4]. The repulsive core is essential not 
only for describing the NN scattering data, but 
also for the stability and saturation of atomic nu- 
clei, for determining the maximum mass of neu- 
tron stars, and for igniting Type II supernova 
explosions [5J. Although the origin of the repulsive 
core must be related to the quark-gluon structure 
of nucleons, it remains one of the most fundamen- 
tal problems in nuclear physics for a long time [6]. 
It is a great challenge for us to derive the nuclear 
potential including the repulsive core from lattice 
QCD. 
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Figure 1. Three examples of the modern NN po- 
tential in 1 Sq (spin singlet and S- wave) channel: 
Bonn[7], Reid93[8] and AV18[9]. 



2. Potentials from Lattice QCD 

2.1. Strategy 

Let me explain the strategy proposed in 
Rcf. |TD] in order to extract the NN potential from 
lattice QCD. 

One first defines a wave function with energy E 
through the equal-time Bethe-Salpeter amplitude 
as 

<p E (f) = (0\N(x,0)N(y,0)\2N,E) (1) 
where r = x — y , \2N, E) is the 2 nucleon eigen- 
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Figure 2. The central potential for the singlet 
NN (red circles) and the effective central poten- 
tial for the triplet NN (blue triangles) as a func- 
tion of r (fm) at ~ 530 MeV in quenched 
QCD 10 . The solid line represents the one pion 
exchange potential (OPEP) given in EafTOl 



state of QCD with energy E, and |0) is the vac- 
uum state. More precisely, the 2 nucleon state 
means a baryon number 2 (or quark number 6) 
state in QCD. We here use a simple local nucleon 
field operator N(x) given by 



N[(x) 



^abc 



(u a (x)C l5 d b (x))q^(x) (2) 



where a, 6, c are color, a is the spinor and / is the 
flavor indices for quark field q(x) = (u(x), d(x)), 
and C = 7274 is the charge conjugation matrix in 
the spinor space. For large r— = \r\, it is shown[?] 
that this wave function behaves as 



A 

kr 



sin(fcr - -1+ 8i(k)) 



(3) 



where ipi is the Z-th partial wave of <p, Si(k) is 
the phase shift of the Z-th partial wave, and k 2 = 
2fiE — iunE with the reduced (nucleon) mass fi 

K). 

Using the above wave function, the correspond- 
ing non-local potential U (x, y) is defined though 
the Schrodinger equation as 



[E-H ]ip E (x) 



d 3 yU(x,y)ip E (y) 



(4) 



where H is the free Schrodinger operator given 
by 

-V 2 



H 



2(i. 



(5) 



Note that U (x, y does not depend on energy E. 

Finally one expands U(x, y) in terms of the 
derivative V as 



U(x,y) = V(x,V)S^(x-y). 



(G) 



Explicitly the leading contribution can be writ- 
ten[?] as 

V{x,V) = V (r) + V a (r)(a 1 -a 2 ) 
+ V T (r)S 12 + V LS {r)L-S + 0{V 2 ) (7) 

where r = \x\, is the spin-operator of each 
nucleon, L = f x p is the angular momentum op- 
erator, S = + 02 ) is the total spin operator, 
and 



Sl2 — 



:(<?! • x)(a 2 ■ x) - (<?i • <t 2 ) 



(8) 



is the tensor operator. Coefficient functions Vx 
has the following decomposition in the flavor 
space: 

Vx = V^) + Vj(r)(ri .f 2 ) (9) 

for X — 0, <7, T, LS, where 7*1,2 is the Pauli matrix 
in the flavor space of each nucleon. 

2.2. First Result 

Applying the strategy mentioned in the previ- 
ous subsection to the quenched lattice QCD cal- 
culation, the first result of the iV./V potential has 
been obtained [TO], as shown in Fig(2l which repro- 
duces qualitative features of NN potentials such 
as the strong repulsive core at short distance ( 
r < 0.5 fm) surrounded by an attraction well at 
medium and long distances. The solid line in the 
figure represent the one pion exchange potential 
given by 



Vv? PEP 



91n (n ■ 7*2) (0=1 ■ <? 2 ) / m 



An 



2/71 j\[ 



x , (10) 

r 

where ~ 530 MeV and m N ~ 1340 MeV 
are taken from the simulation, while the physical 
value of the irN coupling constant, g 2 N /(4n) ~ 
14.0, is used. 
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2.3. Questions and Some Answers 



VcO^SqJPBC v.s. APBC 



600 | 



> 400 



200 



APBC 
PBC 




Figure 3. The central NN potential for the 1 >S'o 
state at E ~ 50 MeV (red bars) and E ~ MeV 
(blue crosses) in quenched QCD. 



So far the strategy to extract NN potentials 
from lattice QCD seems promising. There ex- 
ist, however, a few questions to be considered 
before starting serious investigations of NN or 
more generally baryon-baryon potential using lat- 
tice QCD technique. 

As mentioned before, the potential extracted 
from the wave function is in general energy- 
independent but non-local: 

[E-H ]<peW = J ' d 3 yU(x,y)<p E (y). (11) 

This energy independent non-local potential is 
equivalent to the local but energy dependent po- 
tential, Ve, defined by 



[E-H ](p E {x) 



V e (x)<Pe(x) 



(12) 



The 1st question to our strategy is how Ve{x) 
depends on energy E of the wave function lpe{x). 
The potential in FigJ5J which is obtained at E ~ 
0, can reproduce the correct phase shift at E ~ 
(scattering length) by construction. If Ve(x) were 
strongly dependent on the energy E, it would not 



be so useful: it may not reproduce the correct 
phase shift at E' =/= E. 

The energy dependence of V E has been investi- 
gated in quenched QCD. In FigJ3]the singlet NN 
potentials are compared between E ~ MeV and 
E ~ 50 MeV[?]. For E ~ MeV, the periodic 
boundary condition (PBC) is used, while the anti- 
periodic boundary condition (APBC) for E ~ 50 
MeV. Fluctuations of data with APBC at large 
distances ( r > 1.5 fm ) are mainly caused by 
contaminations from excited states, together with 
statistical errors. A non-trivial part of the poten- 
tials at r < 1.5 fm, on the other hand, are less 
affected by such contaminations. As seen from 
the figure, the NN potentials are almost identi- 
cal between E ~ MeV and E ~ 50 MeV. This 
shows that the potential obtained from our strat- 
egy is useful at low energy at least below 50 MeV 
at this pion mass (m, = 530 MeV) in quenched 
QCD. It will be important to see if this conclusion 
holds for lighter pion mass in full QCD. 

The second question is the operator depen- 
dence of the potential. If the wave function were 
defined in terms of different N(x) other than 
N(x) as 



<p E {f) = {0\N(x,Q)N(y,0)\™,E), 



(13) 



one would obtain the different potential: Ve{t) ^ 
V E (r). Our choice of N(x) correspond to a par- 
ticular choice of the Nucleon effective field theory 
defined by 



-Seff(JV) 



S(N- e abc {u a C lb d b )q c ) . (14) 

Since no field dependence appears at leading or- 
der in the Nucleon effective theory, one expects 
that operator dependence is also weak at low 
energy in QCD. This expectation, however, will 
have to be explicitly checked in numerical simu- 
lations. 

3. Challenges of lattice QCD: Form Quarks 
to Nuclei 

Now I list some challenges of lattice QCD, 
which are relevant in nuclear physics. 
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for determinations of nuclei structures and 
inputs for the equation of state inside Super 
Nova or Neutron Star. 

3. Finally we still need some theoretical under- 
standings for the repulsive core of the NN 
potential. 

In the following I will report some activities 
related to challenges in lattice QCD listed above. 

3.1. Full QCD simulations and hadron 
spectra 



Figure 4. Hadron spectra in the continuum limit, 
together with experimental values. To fix free 
parameters of QCD, m,, vtlk and m~ are used. 



1. Full QCD calculations at physical quark 
masses are definitely important and one of 
the most difficult challenges in lattice QCD. 
Once full QCD configurations are generated 
at physical quark mass in sufficiently large 
boxes, one can try to investigate several in- 
teresting quantities in nuclear physics. For 
example, one can construct a deuteron on 
the lattice, or can extract realistic NN po- 
tentials from lattice QCD. One can also 
make real predictions for hyperon interac- 
tion from lattice QCD. 

2. Calculations of more complicated quantities 
are also challenges in lattice QCD. Not only 
the central part of potentials but also more 
complicated structures of potentials such as 
the tensor force or the L ■ S force should be 
extracted from lattice QCD. The more ad- 
vanced challenge in this category is to de- 
termine the 3 body force among nucleons, 
which can hardly be determined experimen- 
tally. Although one can extract the 3 body 
force from scattering states, it is also a great 
challenge to construct the bound state of 3 
nucleons such as Triton or Helium 3 on the 
lattice. The lattice QCD may supply inputs 




m [lattice unit] 



Figure 5. The ratio rn^/m u d as a function of 
m u d at a = 0.09 fm and L — 2.9 fm. Red circles 
correspond to m„ = 702, 570, 412, 296, 156 MeV 
from right to left. The vertical dashed line is the 
physical point. 



Fig[4] shows the result of hadron spectra from 
one of the most extensive 2+1 flavor full QCD cal- 
culation in Ref. [11) , where both chiral and contin- 
uum extrapolations have been taken. The mini- 
mum pion mass in this calculation is m" lm = 190 
MeV, and m % L > 4 is always satisfied for the 
spatial extension L. The overall scale, the light 
quark mass and the strange quark mass are fixed 
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by Wis, m, and m^. The agreement between lat- 
tice QCD and experiment is excellent, as can be 
seen from the figure. 

It is also reported that the pion mass can be 
smaller than 160 MeV in the 2+1 flavor full QCD 
simulations at one lattice spacing a = 0.09 fm[12j. 
Fig.?? shows the pion mass squired divided by the 
ud quark mass, m^/m u d, as a function of m u d, 
where black solid circles are previous results from 
CP-PACS/JLQCD collaborations [13] obtained at 
heavier quark masses, while red solid circles are 
new ones from Ref.[12j. As quark mass decreases, 
this ratio increasingly deviates from the straight 
line obtained from the fit with black circles. This 
behavior suggests an existence of the chiral loga- 
rithm predicted by the chiral perturbation theory. 
Note however that m™ m L — 2.3 in this calcula- 
tion, where the finite size effect might be sizable 
at the lightest quark mass. To reduce the mag- 
nitude of the finite size effect less than a % level, 
simulations on larger volume, L — 5.8 fm, with 
m„ ~ 140 MeV (m^L ~ 4.1 ), are begun by this 
group. This simulation indeed will realize the real 
QCD calculation. 




Figure 6. 2+1 flavor full QCD results of (effec- 
tive) central potentials at = 702 MeV. 



3.2. Nucleon potential from full QCD sim- 
ulations 

Using PACS-CS 2+1 flavor full QCD config- 
urations, the NN potentials are currently be- 
ing calculated. An example is given in Fig.??, 
where thesingle central potential Vcir; 1 So) and 
the triplet effective central potential V§ s (r; 3 Si) 
are plotted as a function of r at m v — 702 
MeV 14 . By comparing this figure with Fig.??, 
one notices that much stronger repulsive core 
than the quenced case is observed for the poten- 
tials in full QCD. Reason for this difference is 
currently under investigation. 

3.3. Tensor potential 




Figure 7. The tensor potential Vr(r) (red bars) in 
quenched QCD at m T ~ 530 MeV, together with 
the central potentail Vc (r) (blue crosses) and the 
effective central potential V<S ff (r) (green stars). 



Since the spin triplet state with J = 1 contains 
both L — and L — 2, the effective Schrodinger 
equation becomes 



[E - H Q ] <p E (r) 



[V c (r) + V T (r)S 12 }cp E (m 



where S12 is the tensor operator defined before 
and Vr(r) is the tensor potential. Using P, the 
projection operator to 3 5i ( L = ) state, and 
Q = 1— P, the projection operator to 3 Z?i ( L = 2 
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state, the above equation leads to 



PtpE{r) PSi2tpE(r) 
QtfE{r) QSi2<pE{r) 



V c {r) 
V T (r) 



Qy E {r) 

so that the solution can easily be obtained as 

Pip E (r) PS 12 fE(r) 
QfE(r) QS 12 ifiE(r) 

x (E-H )( ^ 



(16) 



V c (r) 
V T (r) 



(17) 



Note that the L-S potential is neglected here since 
it is higher order in the derivative expansion. 
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estimation to the tensor potential in Ref. [?]. In 



Fig. 



dependence of the tensor potential is 



shown. As seen form the figure, the attraction 
of the tensor force becomes stronger as the pion 
mass decreases. One may expect that the tensor 
force becomes strong enough at = 135 MeV 
to have a bound state, deuteron. 

3.4. Hyperon-Nucleon interactions 
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Figure 9. The effective central potential for Ap in 
the 1 5o channel at m % ~ 465 MeV (red circles) 
and Tfi-x ~ 514 MeV (green squares). The inset 
shows its enlargement. 



Figure 8. The tensor potential at m, = 380 MeV 
(green), 529 MeV (blue) and 731 MeV (red). 



In Fig[7J the central Vc(r) and the tensor Vr(r) 
potentials for the spin triplet state, obtained in 
quenched QCD at — 530 MeV, are plotted as 
a function of r, together with the effective central 
potential V^ s (r). One first notice that Vc(r) — 
Vp ff (r). This means that the neglect of the tensor 
potential in the previous calculation [TO] does not 
lead to large systematic errors in the determina- 
tion of the potential. It is interesting to observe 
that the tensor potential seems have no repulsive 
core, and this property is consistence with the 



An important application of the method to ex- 
tract potential from lattice QCD is to investi- 
gate interactions between hyperon-nucleon(YN) 
or hyperon-hyperon(YY), where hyperon(Y) is a 
baryon which has one or more strange quarks. 
These interactions, which are the bases to explore 
the strange nuclear system, must be determined 
experimentally and/or theoretically. Recent sys- 
tematic study for light A hypernuclei( ^H, ^H, 
^He and ^ He), for example, suggests that the 
AN interaction in the 1 S'o channel is more at- 
tractive than the in the 3 Si channel. The present 
YN and YY interactions, however, have large un- 
certainties, since the scattering experiments are 
either difficult or impossible due to the short life- 
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time of hyperons. Therefore it is important to 
determine the YN and YY interactions theoret- 
ically from lattice QCD. 
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Figure 10. The effective central potential for Ap 
in the 3 Si channel at m n ~ 465 MeV (blue tri- 
angles) and ttItt — 514 MeV (magenta diamonds). 
The inset shows its enlargement. 



The NA effective central potentials in quenched 
QCD are shown in Fig[9]for the 1 So state and in 
FigdU] for the 3 Si state at = 465 MeV and 
514 MeV[16j. Note that a clear spin dependence 
of the potential is observed while the pion mass 
dependence is less significant. 



4. Conclusion 

The challenges of lattice QCD, from quarks to 
nuclei, have just begun and will continue for a 
while. New ideas are important, and close col- 
laboration between particle physics and nuclear 
physics is a key for successes of the challenges. 
For example, lattice QCD approach to the po- 
tential should be complemented by more analytic 
methods such as nucleon effective theories [?]. 
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